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, We study holographic superconductors in 3+1 dimensions away from the probe limit, i.e. taking back- 

reaction of the space-time into account. We consider the case of pure Einstein - and Gauss-Bonnet gravity, 
respectively. Similar to the probe limit we observe that the critical temperature at which condensation 
sets in decreases with increasing Gauss-Bonnet coupling. The decrease is however stronger when taking 
back-reaction of the space-time into account. We observe that the critical temperature becomes very small, 
but stays positive for all values of the Gauss-Bonnet coupling no matter how strong the back-reaction of 
the space-time is. 
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The gravity-gauge theory duality pQ has attracted a lot of attention in the past years. The most famous example 
is the AdS/CFT correspondence [2] which states that a gravity theory in a d-dimensional Anti-de Sitter (AdS) 
£Q | space-time is equivalent to a Conformal Field Theory (CFT) on the (d — l)-dimensional boundary of AdS. 

Recently, this theory has been used to describe so-called holographic superconductors with the help of black 
holes in higher dimensional space-time [3J HI [5] and many aspects have been discussed such as holographic su- 
perconductors in Horava-Lifshitz gravity [6 and in Born-Infeld electrodynamics [7] , fermions [8J , the behaviour 
of holographic superconductors in external magnetic fields [5] and at zero temperature [10] , hydrodynamical as- 
pects of holographic superconductors as well as rotating superconductors [12] . Holographic superconductors 
^ ' in extended models that allow for a first order phase transition [T3] as well as holographic superconductors in 
M-Theory [2] have also been studied. Non-abclian (or p-wave) holographic superconductors have been studied 
in[T3[TBl[T71[T51[T51|5niinil21l23]- In QUI H3] a string theory realization of p-wave holographic superconductors 
in the probe limit has been discussed and the Meissner effect has been studied in detail [20]. (For a related 
analytical study see [21].) Sound modes for p-wave superconductors have been considered in [22], while fermions 
in these superconductors have been discussed in [33]. Various other aspects have also been studied [2"4"] . 

The general idea behind holographic superconductors comes from the observation that below a critical 
temperature electrically charged black holes become unstable to form scalar hair, i.e. they possess non-vanishing 
scalar fields on the horizon [3J . The reason for this is that close to the horizon of the black hole the effective mass 
of the scalar field can become negative with masses below the Breitenlohner-Freedman bound [25] such that 
the scalar field becomes unstable and possesses a non-vanishing value on and close to the horizon of the black 
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hole. The value of the scalar field on the AdS boundary is then associated with the corresponding condensate 
in the dual theory. 

In most cases, holographic superconductors have been studied in the "probe limit" neglecting back-reaction 
of the space-time. This limit corresponds to letting the electric charge e tend to infinity or equivalently Newton's 
constant G tend to zero. Backreaction of the space-time was considered in [4] for (2+l)-dimensional holographic 
superconductor. It was found that the qualitative results are similar for small charges, but that suprisingly the 
scalar field can even form a condensate when being uncharged. 

In [26] (3+l)-dimensional superconductors were studied by investigating scalar hair formation on black 
holes in Gauss-Bonnet gravity. This has been extended to higher dimensions in [27] . The motivation for this is 
the apparent contradiction between the Mermin- Wagner theorem that forbids spontaneous symmetry breaking 
in 2+1 dimensions at finite temperatures and the fact that (2+l)-dimensional holographic superconductors 
do exist. Consequently, it has been suggested that higher curvature corrections should suppress condensation, 
where higher curvature corretions can of course only been studied for (3+l)-dimensional superconductors (or 
higher dimensional ones). [261 |2"T] were concerned with the "probe limit" and it was found that condensation 
cannot be suppressed in Gauss-Bonnet gravity. 

In this paper, we are interested in the model studied in |26j but away from the probe limit, i.e. taking back- 
reaction of the space-time into account. We study (3+l)-dimensional superconductors in pure Einstein and 
Gauss-Bonnet gravity, respectively. While for large temperatures, i.e. when the scalar field vanishes identically 
analytic solutions to the equations of motion are known, this is different for a black hole with scalar hair that 
forms below the condensation temperature. These solutions have to be constructed numerically. 

In Section 2, we present the model, the equations of motion and the boundary conditions. In Section 3, we 
discuss our numerical results, while Section 4 contains our conclusions. 



2 The Model 

In this paper, we are studying the formation of scalar hair on an electrically charged black hole in (4 + 1) 
dimensional Anti-de Sitter space-time. The action reads : 

S = J d 5 x^y~ (R-2A+J (R^PR^xp - iR^R^u + R 2 ) + 167rG£ matt er) , (1) 

where A = — 6/L 2 is the cosmological constant and a the Gauss-Bonnet coupling. £ ma tter denotes the matter 
Lagrangian : 

Cotter = -\f mn F mn -(D M ip)* D M ^-m 2 r^ , M,N = 0, 1, 2, 3, 4 (2) 

where Fmn — d^iA^ — dNAM is the field strength tensor and DmiP — diaip — ieAMip is t ne covariant derivative, 
e and m 2 denote the electric charge and mass of the scalar field tp, respectively. 
The Ansatz for the metric reads : 

ds 2 = -f(r)a 2 (r)dt 2 + -±-dr 2 + ^dS 2 (3) 
J{r) L z 



where / and a are functions of r only. The 3-dimensional metric is 

dT,l 3 = { dx 2 + dy 2 + dz 2 for k = (4) 



dE 2 for k = — 1 



where k denotes the curvature of the 3-dimensional space. We are only interested in plane-symmetric black 
holes in this paper, so we will set k — 0. However, we will keep the k in the equations for completeness. 
For the electromagnetic field and the scalar field we choose [4] : 

A M dx M = <t>(r)dt , ^ = V(r) (5) 
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such that the black hole possesses only electric charge. 

The coupled Einstein and Euler-Lagrange equations are obtained from the variation of the action with 
respect to the matter and metric fields, respectively. They read : 
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(7) 
(8) 
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where 7 = lQirG. Here and in the following the prime denotes the derivative with respect to r. In [26 , 
(3 + l)-dimensional holographic superconductors have been studied in the probe limit corresponding to 7 = 0. 
For 7 / we take back-reactions of the space-time into account. Note that this limit is equivalent to letting 
e — > 00 since we can preform the rescalings %p — > ip /e, <j> — » <p/e and 7 — > e 2 7. Hence without loosing generality 
we can set e = 1. 

In order to find an explicit solution of the equations of motion, we have to fix appropriate boundary condi- 
tions. In the following, we are interested in the formation of scalar hair on electrically charged black holes with 
horizon at r = ru such that 

f(r h ) = (10) 
with a{rh) finite. In order for the matter fields to be regular at the horizon we need to impose: 



(r h ) = , V'M 



m 2 il) (r 2 + 2ak) 



2rk + Ar/L 2 - 7 r 3 (m 2 ^ 2 + c/)' 2 /(2a 2 )) 



(11) 



Asymptotically, we want the space-time to be that of a Reissner-Nordstrom-Anti de Sitter black hole, i.e. we 
can choose a(r — > 00) — > 1. Other choices of the asymptotic value of a(r) would simply correspond to a rescaling 
of the time coordinate. The matter fields on the other hand obey [26] : 



(r > 1) = n - p/r 2 



^ r » !) = — + ^ 



(12) 



with 



A- = 2 — \/4 — 3(L c ff / L) 2 , A+ -2 + v /4-3(L cff /L) 5 



T 2 



2q 



1 - y/l -4a/L 2 



^ L 2 (l-a/L 2 + 0(a 2 )) . 

(13) 

Note that the value of the Gauss-Bonnet coupling a is bounded from above : a < L 2 /4 where a = L 2 /4 is the 
Chern-Simons limit. For larger values of a the solution would possess a naked singularity. 

The parameters p, p are the chemical potential and density of electric charge, respectively. Along with [2 6) 
we choose ip- = 0. will correspond to the expectation value (O) of the operator O which in the context of 
the gauge theory-gravity duality is dual to the scalar field and as such represents the value of the condensate. 

There are analytic solutions of the equations of motion for ij)(r) = : 



f(r) 



k + TT 1- 
2a \ 



4a 

l -T 2 



AaM Aajp 2 



a{r) = 1 



(0 = 4-4 



(14) 



where M and p are arbitrary integrations constants that can be interpreted as the mass and the charge density 
of the solution, respectively. In the limit a — > 0, the metric function f(r) becomes f(r) = k + jj — + 
These solutions are electrically charged black holes which are the only solutions for temperatures larger than 
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the critical temperature T c . For T < T c these solutions will be unstable to form scalar hair, i.e. develop a 
non-vanishing value of tp on the horizon. The aim of this paper is to study the formation of scalar hair black 
holes in dependence on a and 7. The temperature mentioned here corresponds to the Hawking temperature of 
the black hole and reads 

= ^-f(r h )a(r h ) , M, N = 1, 2, 3, 4 . (15) 

47T 

In the gauge theory - gravity duality T c is the temperature below which superconductivity appears. 



T= - r V-9 tt 9 MN dMgttd N gu 



3 Numerical results 

In the following we are only interested in the plane-symmetric black holes with k = 0. The equations of motion 
©-© depend in principle on a number of constants but due to the scale invariances noted in [3] two of them 
can be scaled out and hence be fixed to particular values without loosing generality. In the following we fix 
77j = 0.5 and L = 1. Along with [26] we set m 2 = — 3/L 2 = — 3 which guarantees the stability of AdS^ since 
m 2 < fTigg = — 4/L 2 with m| F the Breitenlohner-Freedam mass [2"5] . 

To find a unique solution to the equations of motions we fix the boundary conditions (|10[) . (jlip at the 
horizon, and choose a(oo) = 1, tp- = 0. In addition we fix -0+ to a particular value. In this way we are able 
to construct branches of solutions labelled by the parameter Note that p and p are uniquely fixed by 
the choice of ip+ and are not free parameters. However, in the literature on holographic superconductors, the 
solutions are typically presented for fixed electric charge density p, while the horizon value is treated as a 
free parameter. These two approaches are connected to each other. Indeed, it is easy to convert a branch of 
solutions with fixed and varying ip + - where p = p{ip+) - into a branch of solutions with constant charge 
density. In the following we denote quantities corresponding to a fixed charge density by a hat. Setting p = 1, 
the relevant Hawking temperature T and condensate ip+ & re respectively given by 

f = ^N ' ^ + = J^m > (16) 

where A+ is defined in (| 13)) . 

3.1 Effect of back— reaction in Einstein gravity 

This corresponds to the case 7 7^ and a = 0. We solved the equations for several values of 7 and -0+ and find 
that solutions exist for generic values of these parameters. 

When studying solutions for 7 fixed and varying ip + we find that in the limit 0+ — > the solutions tend to 
(|14|) for very specific values of p and p which depend on the choice of 7 and can only be determined numerically. 
Correspondingly, the critical temperature T c at which -0 + = can also only be determined numerically. We 
find the values given in the table below : 
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0.0 


0.025 


0.05 


0.1 


0.15 


0.2 


0.3 


0.35 


47rT c 


2.48 


2.02 


1.61 


0.99 


0.57 


0.33 


0.10 


0.06 



For large 7 the construction of the solutions becomes increasingly difficult. In principle we would want to 
know what happens for very large 7. In order to understand this, we fitted the numerical data and found that 

T c w 0.198 -exp (-10.6- 7) p 1/3 (17) 

fits the data for 7 > 0.2 very well. This on the other hand means that no matter how large we choose 7, we 
will always have T c > 0. This has already been observed for superconductors in (2 + 1) dimensions [4], where 
it was shown that T c > in the limit e — > which corresponds to 7 — > 00 here. Apparently, this phenomenon 
persists for (3 + l)-dimensional superconductors. 
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Fixing 7 and increasing the value of the condensate ip + we find that the values a(r/j) and ^'(r^) slowly 
approach zero. At the same time, the function /(r) develops a local maximum and a local minimum at values 
Hvi, fm such that rh < j~m < r m < oo- This is illustrated for the metric functions f(r) and a(r) in Fig. [TJ (left) 
for 7 — 0.2 and three different values of 
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Figure 1: (left) The metric functions f(r) and a(r) (left) and the matter functions cj>(r), ip(r) and 4>'(r) (right) for 
7 = 0.2 and for three values of ip+. 

This would suggest that for sufficiently large -0+ the function f(r) develops a double zero at r = r m which 
would correspond to the formation of an extremal black hole with vanishing Hawking temperature. A detailed 
analysis however shows that the value /(r m ) remains strictly positive, while the value of a(r/j) decreases with 

increasing according to an exponential behaviour a(rh) ~ exp (— cip+) with c some constant. This result 
suggests that the black hole solutions are not limited by a maximal value of the condensate an d that the 
temperature stays positive for all values of 0+. 

In FigfT] (right) we show the matter field functions <j>{r), ip(r) and 4>'{r). We observe that for fixed 7 and 
increasing ip + the maximum of 4>'(r) is pushed further away from the horizon of the black hole. Indeed, for 
small values of 7 and V'+j the maximum of 4>'(r) is on the horizon of the black hole. Note that when fixing ip + 
and increasing 7 we observe a similar phenomenon. 

To understand how the scalar field ip becomes unstable close to the horizon, we plot the effective mass 



2 A 2 U 

e A t g = - 



3 

1? 



fa 2 



(18) 



in Fig 12] for L — 1. Indeed, the effective mass drops below the Breitenlohner-Freedman bound of — 4/L 2 = —4 

close to the horizon. For fixed ip+ and increasing 7 the quantity —j? — -^2 becomes more narrow and smaller 
in absolute value. For fixed 7 and increasing ip + it becomes broader and larger in absolute value. 

To understand this in more detail note that close to the horizon the functions can be expanded as follows : 
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Figure 2: The value of the effective mass m 2 cfi = m 2 + e 2 A 2 g u 
several values of 7 and tp+. 
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close to the horizon of the black hole for 



Moreover note that there are the following relations between the values of the functions at r — : 



4>"{r h ) = 0'(r h ) 
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f'(r h ) = 4^ - r hl |mV(r k ) 2 + 



^'(r h )f'(r h ) = m 2 il;(r h 



L 2 



2a(r h y 
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a'(r h ) = r h -f [ a(r h )tp' '(r h ) 2 + e 



.<f>'(r h ) 2 i>(r h Y 



(23) 

(24) 
(25) 

(26) 



a{r h )f>{r h y 

Several quantities characterizing the solutions are given in Fig|3] for 7 varying and ip+ — 0.2 fixed (left) 
and for tf> + varying and 7 = 0.2 fixed (right), respectively. First note that ip(rh) increases with 7 and tp+, 
respectively. That ip{ r h) is an increasing function of the condensate "0+ was already noticed in |26j for the 
probe limit. Here, we find in addition that the stronger the back-reaction the higher the value tp{rh) for a given 
condensate ?/>+. Since we would like -0 to have its maximal positive value on the horizon we have ip{ r h) > 0, 
"0'( r /t) < and then from (|24]l obviously f'(rh) > 0. For 7 = 0, the value of f'(rh) = ^jh which for = 0.5 and 
L = 1 is just f'(rh) — 2.0. For increasing 7 the value f'(rh) is first decreasing due to the decrease of the electric 
field 4>'(rh)/a(rh) on the horizon (the negative term in (|25[1 ) and then for sufficiently strong back-reaction ip( r h) 
becomes larger and larger such that /'(r/,) starts increasing again. This is similar for 7 fixed and varying ip+. 
For ip + — 0, the solution is given by (fT4|) . For ip + > but small the electric field on the horizon 4>'(rh)/a(rh) 
decreases and leads to a slight decrease in /'(r/j). For increasing ip + the value ip( r h) becomes larger and f'(rh) 
increases. For both 7 and ip+ increasing, respectively, the value of a(rh) decreases from a(r/j) = 1. 

Finally, the charge density p decreases for small ip + and increases for larger ip + when 7 is fixed. For fixed 

and varying 7 the behaviour is qualitatively similar. 

When studying the dependence of the condensate ip+ on the temperature T one can take two different 
viewpoints. First we consider the system for fixed horizon (r^ = 0.5 here). The dependence of < O >= tp+ as 
a function of the temperature T is given in Fig. |4] (left) for several values of 7. 
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Figure 3: Several quantities characterizing the black holes in pure Einstein gravity (a = 0) with tp+ = 0.2 for 7 varying 
(left) and for 7 = 0.2 and ip+ varying (right). 

The solid line represents the condensate ip+ and the dashed line represents the charge density p, respectively. 
For completeness we give the corresponding lines for 7 = (dashed line) , where the bullet indicates the minimal 
value of p. Apparently, the behaviour is quite different when comparing large and small 7. For small 7 the 
temperature of the condensate is larger than the temperature of the critical limit, for large 7 it is vice versa. 

Following the literature, we also present our results for fixed charge density p = 1. The dimensionless 
quantity < 6 > 1/A + /f c = (V3+) 1/A +/T c (with A+ = 3 in the limit a = 0) as a function of the rescaled 
temperature T/T c is given in Fig. |4] (right). Qualitatively, the behaviour for large 7 is similar to that for small 
7. However, the condensate can become quite large when increasing 7. Moreover, the critical temperature T c 
at which = decreases with increasing 7. 

Note that we are in fact plotting (e^'+) 1 / A +/T c but that e doesn't appear here due to our choice e = 1. 
Moreover, comparing our results to the (2+l)-dimensional case [4] our choices of 7 = 0.0, 0.05, 0.1, 0.2 and 0.3, 
respectively would correspond to e = 00, 4.47, 3.16, 2.24 and 1.83 when setting 7 = 1 instead of e = 1. 

3.2 Effect of back— reaction in Gauss-Bonnet gravity 

This corresponds to the case 7^0 and < a < 0.25. The 7 = limit was studied in [26]. Our numerical 
results indicate that also for 7 ^ the presence of the Gauss-Bonnet term leads to a decrease in the critical 
temperature. This is shown in Fig. [S] where we give the dependence of T c on a for different values of 7. 
Apparently also back-reaction on the space-time cannot suppress condensation, i.e. the critical temperature 
stays positive for all values of 7 and a that we have studied in this paper. We find e.g. for a — 0.1 that 
T c = 0.185/9 1 / 3 for 7 = (in agreement with [26]), T c = 0.051P 1 / 3 for 7 = 0.1 and T c = O.OO8P 1 / 3 for 7 = 0.2. 
For a = 0.25 we find T c = 0.158P 1 / 3 for 7 = (again in agreement with [26]), T c = 0.024/5 1 / 3 for 7 = 0.1 and 
T c = O.OOlp 1 / 3 for 7 = 0.2. Hence, the critical temperature can become arbitrarily close to zero for a and 7 
large enough, however within our numerical accuracy, we never find T c = for finite values of 7 and a < 0.25. 

4 Conclusions 

In this paper, we have studied holographic superconductors in 3+1 dimensions away from the probe limit. We 
considered the case of pure Einstein and Gauss-Bonnet gravity, respectively and have constructed numerically 
electrically charged black holes that carry scalar hair. For pure Einstein gravity we find that in agreement with 
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A A 

T T/T c 

Figure 4: The values and p as functions of the temperature T for several values of 7 with rn = 0.5 (left). The 
values (t/>+) 1/ ' 3 /T , c as function of T/T c for several values of 7 with p = 1 (right). 

the results for holographic superconductors in 2+1 dimensions [4] the critical temperature at which condensation 
sets in is strictly positive for all values of the gravitational coupling. Considering Gauss-Bonnet corrections 
further decreases the critical temperature, but all our numerical results indicate that it stays positive when taking 
back-reaction into account. Hence, even when taking the gravitational coupling to infinity - which corresponds 
to letting the electric charge e of the condensate tend to zero - there would still be condensation. Similar to 
2+1 dimensions this signals the existence of an additional instability of the scalar field. The explanation is 
similar to that in 2+1 dimensions [4j [5] : since the scalar field is uncharged the instability cannot be caused by 
the spontaneous symmetry breaking. Rather it is caused by the fact that for large 7 the black hole is close to 
the extremal limit in which its horizon geometry would correspond to AdS2 x K3. In AdS2 the Breitenlohner- 
Freedman bound [25] is m 2 BF = — 1/(4L 2 ). Hence a scalar field with mass m 2 — —3/L 2 that is stable in AdSs 
is certainly unstable in AdS2- 

When considering non-abelian holographic superconductors it has been observed that the phase transition 
that leads to the formation of vector hair becomes first order if the gravitational coupling is large enough [T7] . 
It would be interesting to see how the Gauss-Bonnet term influences this result. 

Acknowledgments We thank J. Erdmenger for comments on this manuscript. YB thanks the Belgian FNRS 
for financial support. 
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